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ABSTRACT 



The direction of sand movement along the Carmel River State Beach 
was qualitatively determined by diving observations, a bathymetric survey, 
wave refraction diagrams and a sediment size analysis of 18 samples. 

The primary source of sediments for the beach appears to be the 
Carmel River which flows only seasonally. Sedimentary material is intro- 
duced into the bay after winter precipitation provides a sufficient 
amount of run-off to warrant the opening of the river mouth by bulldozer. 

The fine sedimentary material is lost offshore and the coarser ma- 
terial is either redeposited on the beach or is carried south with the 
littoral drift and deposited at the nodal point in the sand transport 
pattern. This node is located on the northern edge of the head of the 
Carmel Submarine Canyon. 

Winter storms probably induce slumping or gravity sliding and much 
of the material is carried to deeper water by the canyon. 
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I . INTRODUCTION 



A. OBJECTIVE 

This study was undertaken in order to determine the direction of 
sand movement along the Carmel River State Beach, Carmel, California, 
and provide an insight into the methods by which this sand is introduced 
into the Carmel Submarine Canyon. In order to accomplish these objec- 
tives, the area of study was examined by the author using SCUBA, a bathy- 
metric survey was conducted, wave refraction diagrams were constructed, 
an attempt was made to measure bottom currents, and sediment samples were 
taken for textural analysis. 

B . AREA DESCRIPTION 

Carmel Bay, located about 9 km south of Monterey Bay on the central 
California coast (Fig. 1 ), is an embayment that contains the head of 
the Carmel Submarine Canyon (Fig. 2). It is a very open, rectangular- 
shaped bay with dimensions of k-3 km by 3.6 km in the north-south direc- 
tion and east-west direction, respectively. 

Neighboring communities include Pebble Beach and Carmel, California. 
At the southern end of the bay is Point Lobos State Reserve. Two water- 
courses, the Carmel River and San Jose Creek, drain a land area of 
approximately 6?0 km^ and terminate near the head of the Carmel Canyon. 
[California State Department of Water Resources, 1969]* 

The study area (Fig. 3 and lj) is bounded by Carmel (Abalone) Point 
to the north and San Jose Creek to the south. A depth of 73 ft was 
arbitrarily chosen for the western boundary and the beach face formed 
the eastern boundary. 
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Figure 2 . 



Carmel Bay Bathymetry (Contour Interval 
(After Zardeskas, 1 971 ) . 
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Figure li. 



Aerial Photograph of Southern End of 
Carmel Bay Showing the Area of Study. 
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Figure 5- Bulldozer Opening the Mouth of the Carmel River- 
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The primary source of sediments for the study area was assumed to be 
the Carmel River. Sand is introduced by the river primarily during the 
winter and spring when rainfall is sufficient to maintain a flow. During 
the summer and fall, the beach closes the river and a brackish lagoon 
forms behind the beach, thus ceasing the river-carried sand supply. A 
secondary source of sediment is the shoreline -outcropping rock. 

When the winter rains have supplied a sufficient amount of precipi- 
tation to produce a steady flow in the river and the lagoon is filling, 
the mouth of the river is opened by a bulldozer (Fig. 5>)* Failure to 
open the river at the proper time could result in flooding of homes at 
the lower end of the Carmel Valley. This operation is performed under 
the direction of the head ranger at Point Lobos State Reserve and normal- 
ly occurs once or twice each winter. During the winter of 1 966-1 967 , 
this operation was repeated 25> times because of alternating rainy and 
dry periods. The possibility exists that heavy rains coupled with heavy 
surf, a high tide, and a strong onshore wind might prevent this opening 
and flooding would occur [Mr. Don Rich, personal communication] . 

One cut is all that is required as the flow of the river widens the 
cut until the river is at least 100 ft wide at the beach (Fig. 6). The 
author estimated that, at the point of maximum opening, 10,000 to 20,000 
cubic yard of sediment had been removed from the beach. 

Much of the fine sediment material was observed moving out into the 
bay and was obviously lost to the study area. The coarse sediment un- 
doubtedly settled quickly and was either moved with the littoral drift 
or redeposited on the beach by wave action when the river mouth was 
closed again. 
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Figure 6. Mouth of the Carmel River After Flow Has Widened 
the Bulldozer Cut. 
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C . PREVIOUS WORK 



The geology of the area in the vicinity of Carmel Bay has been de- 
scribed in some detail by Lawson [ 1 893 ] ^ Beal [l9l£], Bowen [l96£], 
Nili-Esfahani [l 963>] and Simpson [l 972] . Of interest to this study are 
the formations that are in the study area or are associated with the 
Carmel River and San Jose Creek waterbeds . The Santa Lucia porphyritic 
biotite granodiorite (Cretaceous) is apparently the most important source 
of sediments for the area. The Carmelo Series (Paleocene), interbedded 
pebble conglomerate, sandstone, and siltstone, and Monterey shale (Upper 
Miocene), a silicious shale, provide a smaller amount of sediments. Sur 
Series metamorphic gneisses (pre-Cretaceous) apparently contribute sig- 
nificantly to the heavy minerals supplied to the area as they dominate 
the upper regions of the Carmel River watershed [ California State Depart- 
ment of Water Resources, 1 96 9 ] - 

Diver description of the study area has been almost nonexistent. 
Numerous dives have been conducted in and around the head of the Carmel 
Submarine Canyon. Shepard and Dill [1966], McLean and Peckham [1961], 
Moritz [1968], and Wallin [1968] have provided descriptions of San Jose 
Creek beach and the head of the canyon. Sand falls and slump scars have 
been noted by Moritz . 

Bathymetric surveys of the bay were accomplished by the U.S. Coast 
and Geodetic Survey (USCGS) in 1 883 (No. l5U8a) and 1933 (No. H9U93)- 
Shepard and Emery [l9hl] have conducted a lead line survey of the head 
of the Carmel Submarine Canyon. Zardeskas [ 1 971 ] conducted the first 
bathymetric survey using electronic navigation and echo-sounding equip- 
ment. This survey is considered to be more accurate than previous sur- 
veys but docs not provide sufficient detail inside the 1 O-fm contour. 
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Studies of the currents in the bay are almost non-existent. Tides 
of the bay are the mixed type that is characteristic of the Pacific Coast. 
Diurnal differences between mean higher high water and mean lower low 
water average 5-2 ft [U.S. Coast and Geodetic Survey, 1 971 ] . 

Bascom [196I4] conducted a study of the seasonal changes in the pro- 
file of the Carmel River State Beach. These seasonal changes were of 
little significance in this study. However, these seasonal cycles show 
that sand is moved from the beach to an offshore bar by winter storms 
and is carried back onto the beach by summer swell. 

Carter [ 1 971 ] studied the marine sediments of the bay using grain- 
size analysis techniques and Griffin [l 969] has conducted a heavy-mineral 
analysis of the beaches in the bay. 
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II. METHODS OF DATA ACCUMULATION 



A. BATHYMETRY 
1 . Survey 

A bathymetric survey of the area was conducted on May 12, 1972. 
This survey was undertaken in order to provide information for compari- 
son with previous surveys by the USCGS and the Naval Postgraduate School 
(NPS). Large depth discrepancies between the various surveys could pro- 
vide an indication of slumping. Additionally, the survey provided an 
accurate nearshore chart for wave refraction diagrams. 

Four transit locations along the beach in the study area labelled 
A, B, C, and D (shown on the finished chart. Fig. 13) were located using 
the U.S. Geological Survey Horizontal Control Station "Whaler's Knoll 
Number 2" (Fig. 2), the new Carmel Mission Spire, and the references 
used in the 1933 USCGS survey that were designated "Hudson" and "Granite 
Point" on the finished chart. Two Dietzgen 6000 Series land surveying 
transits were used, first at locations A and B and last at locations C 
and D, to determine the position of the survey boat. A large marker was 
placed on Whaler ' s Knoll Number 2 which could be seen from the transit 
locations. This served as a zero reference for the transit operations 
and all angles were measured relative to this marker. 

An Apelco Model MR-201 B fathometer was mounted on the transom of 
a 1 3 -ft small craft powered by a 9-5~hp outboard motor. An aluminum 
pole painted international orange was mounted vertically above the fath- 
ometer transducer. Communication between the survey boat and the transit 
operators was maintained with 1 -w two-way radios at each location. The 
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position of the survey boat was determined at least every minute during 
the survey. 

Track lines were intended to be perpendicular to the beach at 
30-yard intervals. This was difficult to maintain because of large 
amount of kelp ( Macrocystis pryif era ) and difficulty in seeing prelocated 
beach reference marks. 

2 . Depth Corrections 

All depths on the finished chart were corrected to mean lower 
low water datum. In order to correct for the depth of the transducer, 
differences in calibrated sound velocity, and instrument error, compari- 
sons between the fathometer readings and lead-line readings were made 
periodically over a known flat bottom. 

A tidal correction was applied using the NPS tide gage located 
on Monterey Municipal Wharf Number Two in Monterey Bay (just north of 
Carmel Bay) as a reference. Tidal differences between Monterey Bay and 
Carmel Bay rarely exceed 0.3 ft [Dr. W.C. Thompson, personal communica- 
tion]. Since the fathometer sound beam spreads spherically within its 
beam width, a correction is required if the bottom is not horizontal. 

True slope of the bottom was obtained by plotting the various 
data points and measuring the distance between two points across the 
slope . The fathometer readings are known at each point and the slope is 
easily obtained. A correction for slope was applied using tables prepared 
by Zardeskas [ 1 9 71 ] . 

B . DIVER OBSERVATIONS 

An examination of the study area was made over a 7-month period 
using SCUBA, during which 23 dives were made in the area. A log of six 
of the more interesting dives is contained in Appendix A. On all dives 
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an effort was made to record general bottom types, sand thicknesses, 
unusual features, bottom slopes and, when applicable, sediment and rock 
samples were taken. 

C. LITTORAL DRIFT 

Littoral drift was qualitatively determined using statistical wave 
data from National Marine Consultants (NMC) [i960], and wave refraction 
methods given in U.S. Naval Hydrographic Office [1958] Publication Number 
H.O. 23b' No attempt was made to quantitatively determine the longshore 
sand transport because of insufficient wave data and because empirical 
relationships to determine volume transported are based on a smooth and 
straight beach [Bowen and Inman, 1966]. The beach in the study area is 
broken in numerous places by large granodiorite outcrops. 

1 . Statistical Wave Data 

The wave data used are the results of hindcasting from weather 
map analyses for the three years 1956, 1957 and 1958 (Fig. 7). Sea and 
swell data are given by period, deep water wave height, frequency of 
occurrence and direction. Direction is given in 22.5° increments 
(¥, WNW, NW, etc.). 

None of the seven stations along the California Coast for which 
the data are presented were located in the vicinity of Carmel Bay. The 
data of the nearest station to the north (Number 3) } located approximate- 
ly 50 nm west of the entrance to San Francisco Bay (Latitude 37*6° N, 
Longitude 123-5° V/), and the nearest station to the south (Number U) 
(latitude 35-5° W, Longitude 122.0° W), showed no significant differences. 
For this reason, the data from Station I4 were used. 

The NMC information is incomplete because it is based on synoptic 
meteorological data from the North Pacific and does not consider: 
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Figure ?• Swell Data from National Marine Consultants 
for Station U. (after Cherry, 1969). 
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a. local waves generated by the diurnal sea breezes, and 

b. southerly swell. 

The first, local waves, could be of considerable significance in 
sand movement. The second, because of the great amount of refraction 
that would be involved, would probably not contribute significantly to 
the sand transport. 

The data show that over 80% of the deep water wave energy that 
arrived off the coast comes from the directions of west, west-northwest 
and northwest. Over b0% of the total energy comes from the northwest 
alone [ Cherry, 1 965 ] • 

2 . Construction of Wave Refraction Diagrams 

Wave refraction diagrams were constructed for waves from the 
west, west -northwest and northwest by the author and by previous inves- 
tigators [Wallin, 1968]. The methods given in H.O. 23ii were used. All 
calculations were accomplished using a computer program that eliminates 
the need for tables to determine the ratio of the depth of the water to 
the deep water wave length. A copy of the computer program, formats for 
data cards and a sample of the output are listed in Appendix B. 

Direction of transport for San Jose Creek Beach, located just 
south of the study area, and Carmel River State Beach were determined. 
Refraction diagrams were considered more accurate than previous diagrams 
because of the improved bathymetric data provided by Zardeskas [l97lj 
and the author's own nearshore survey. 

A wave refraction diagram for swell from the west-northwest with 
a period of 11 sec is shown in Fig. 8. 
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Figure 8. Wave Refraction Diagram for West -Northwest Swell with 
a period of 1 1 sec . 
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D. CURRENT STUDIES 



1 . Savonius Rotor Current Meter 

Two attempts were made to measure bottom currents outside the 
surf zone . The first involved the use of a Hydro Products Model 501 in 
situ current recording system. This instrument measures current speed 
with a Savonius rotor and current direction with a direction vane. Di- 
rection, speedy and temperature are recorded with a Rustrak recorder. 
Current speed in shallow water can be measured with this sytem only if 
the water surge component is negligible compared to the current [Dr. E.B. 
Thornton, personal communication]. For this reason, the meter was operated 
for 1 week in the spring with the hope that 1 or 2 days of minimal wave 
activity might ensue. 

The meter required a vertical orientation in order to eliminate 
any gravity component on the directional vane when recording. Since the 
bottom slope was 1 i|° at the chosen location, the vane was fitted with a 
small glass floatation chamber to eliminate this problem when submerged 
in salt water. Additionally, the meter was equipped with a plexiglass 
cover for diver observation of the recorder. 

A 550-lb rectangular (18 inch x 18 inch x 2 inch) lead base im- 
bedded with stainless steel studs for meter attachment was prepositioned 
using the R/V ACANIA in 60 ft of water in the center of the study area. 

One of the sediment samples (Number 6, Fig. 11) was taken next to this 
base. The meter was lowered from a small boat and bolted to the base by 
divers (Fig. 9). 

2 . Ducted Current Meter 

The second effort used a Marine Advisor's ducted current meter 
system attached to a base fabricated at the machine facility of NPS 
(Fig. 10). This system, unlike the Savonius rotor model, can essentially 
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Figure 10. Diver Placing Ducted Current Meter and Base on the Bottom. 



eliminate the wave surge component simply by physically aligning the 
axis of the meter parallel to the crests of the wave-induced sand ripples. 
Its major disadvantage is the requirement for surface recorders. The 
power for the recorders was supplied by a 12-v automobile battery and a 
DC/AC inverter. The recorders consisted of a ducted current meter read- 
out Model S— 6A (0-1 5 knot range) and a Varian Model G-11A strip chart 
recorder adjusted for a full scale deflection of 2 knots. The meter was 
placed in 60 ft of water directly seaward of the Carmel River lagoon. 

The axis of the ducted meter was oriented so that the wave surge component 
was minimal. Data was taken for periods of approximately 30 min.. 

E. SEDIMENT SIZE ANALYSIS 

1 . Collection and Analysis Procedures 

Eighteen sediment samples were taken in the area and grain-size 
analyses were conducted using the procedures outlined in Krumbein and 
Pettijohn [ 1 938 ] . Twelve of the samples were taken by the author using 
SCUBA, five were taken at selected locations along the beach, and one 
(sample Number 18) was taken in the bed of the Carmel River in the vi- 
cinity of the U.S. Geological Service River Gaging Station near the 
bridge of California Highway 1 (Fig. 11). All diver samples, with the 
exception of sample Number 12, were taken across the crests of sand 
ripples to insure uniformity (Fig. 12). Samples Number 11 and 12 were 
taken in a crest and adjacent trough respectively for comparison. 

From each sample a subsample of approximately 1)0 g was taken. 

The subsample was washed with distilled water in order to remove the 
salt, allowed to settle, and the excess water was decanted. The sub- 
sample was then wet sieved through a stainless steel 1|0 screen. No fine 
fraction ( > 1)0) of any significance was noted for any of the samples. 
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Figure 11. Locations of Sediment Samples. 
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Figure 12. 



Diver Taking Sediment Sample. 








The coarse fraction was placed in a beaker , dried and weighed. After 
weighing, the coarse fraction was dry sieved through 8-inch U.S. stan- 
dard screens at 0.50 intervals using a Ro-Tap automatic shaker. The 
fraction retained on each of the screens was then weighed to the nearest 
0.1 mg. Finally, the cumulative weight of fractions was compared to the 
original weight to check for gross errors. 

2 . Computer Analysis 

Folk and Ward [1957] statistical parameters were calculated 
using an IBM 360 computer and a computer program prepared by W. R. 
Anikouchine [Dinger, 1970 , p . 31 ] and slightly modified by Miss Sharon 
Raney (NPS Computer Facility). A copy of the computer program, formats 
for data cards and a sample of the output are listed in Appendix B. 
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III. RESULTS 



A. BATHYMETRY 

Chart preparation was accomplished by first preparing a smooth sheet. 
This required the plotting of all fixes that had been taken on a 1 :J?000 
scale sheet, printing the corrected depths on the plotted positions and 
contouring the sheet in increments of 10 ft. The contours range from 
20 ft to 100 ft. The 10-ft contour was omitted because of insufficient 
data. Moderate surf activity on the day of the survey precluded boat 
operations very near shore. The shoreline for the chart was obtained 
from the USCGS 1933 smooth sheet. 

The finished chart (Fig. 13) was prepared on tracing paper from the 
smooth sheet. It has been reduced for inclusion in this paper. 

Of interest is the shoreward trend of the contour lines between 
transit locations B and C. This might possibly be due to a deprivation 
of sediment supply caused by a large granodiorite outcrop. Additionally, 
there exists a slight seaward extension of contour lines just offshore 
of transit location C. This might be caused by a larger rate of sediment 
supply and/or an underlying formation. This seaward bend of contour 
lines is dramatically shorn in Fig. 2. 

Comparison of the finished chart with the USCGS 1883 and 1933 surveys 
revealed no significant depth changes. The study areas of this chart and 
that of Shepard and Emery [l9hl] did not sufficiently overlap and thus 
provided little useful correlation. The 1971 survey (Fig. 2) did not, 
because of limitations of the towed sounder, survey this close to shore. 
The 10-fathom curve of the 1971 survey was taken from the USCGS 1933 
survey . 
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Figure 13- Nearshore Bathymetric Chart. (10-ft contour omitted because 
of insufficient data). 
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Shepard and Emery had previously compared their survey and the USCGS 
1883 survey and shallower soundings near the canyon head (south of the 
study area) indicated that a slide or a slump had occurred. 

B. DIVER OBSERVATIONS 

The area of study was broken into Subareas 1 through 7 (Fig. 11*) in 
order to describe observations by the author using SCUBA. 

Subarea 1 was entirely granodiorite from the shoreline to a depth of 
65 ft. There appears to be no way that sediments could be transported 
across this barrier. Wave refraction studies [Dr. W.C. Thompson, per- 

t 

sonal communication] and heavy mineral analysis [Griffin, 1969] reinforce 
this assumption. 

Subarea 2 is almost entirely covered with a thick layer (greater than 
5 ft) of fine to coarse sand. One small granodiorite outcrop exists in 
approximately 25 ft of water just west of the Carmel River. Bottom slopes 
in the area range from nearly horizontal in shallow areas to 25 ° -30° in 
the southwestern corner of the area. 

Subarea 3 is a combination of rock and sand bottom. There is no 
physical barrier to sand movement in any direction in this area. 

Subarea lx is entirely sand with the exception of one large granodio- 
rite outcrop. This outcrop is narrow (25 ft to 50 ft) and long (approxi- 
mately 100 yards) and is oriented parallel to the beach. On the western 
side of this unbroken outcrop, bottom slopes range from 5° to 20°. The 
bottom on the eastern side has very gentle slopes. Sand sizes on the 
western side are noticeably finer than on the eastern side. This large 
outcrop precludes any significant sediment transport in a westerly 
direction. Sediment movement in the area appears to be confined to the 
area between the beach and the outcrop and is generally in a southerly di- 
rection. The grain size of the sand gave evidence of much local erosion. 
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Subarea 3 contains the southern end of the granodiorite outcrop de- 
scribed in Subarea 1*. The outcrop is broken in several locations near 
the end. Through these breaks there was evidence of sediment movement 
as large fans or aprons of sedimentary material existed at the base of 
each of the breaks. Some sediment movement was observed by the author. 
It is presumed that a large swell is required before large amounts of 
sediment are moved. The remainder of the area is covered with a very 
thick layer of very coarse sand. 

Subarea 6 is located just seaward of a large rock outcrop that is 
west of the entrance on San Jose Creek. This outcrop is awash about 3 
to U ft at low tide and is called the "washrock" by local divers. This 
area contains the head of the Carmel Canyon. Slopes in this area are no 
less than ij.0 0 . The bottom, although sides would probably be a more 
appropriate term, is entirely made up of large granodiorite boulders. 
Depths in the area reach 120 ft within a few tens of yards from the 
"washrock" . 

Subarea 7 is almost entirely rock covered. There appears to be 
little or no sediment transport across this area. 

Subarea 3 is the most interesting area observed. Previous studies 
have described "sand chutes" along the north wall of the canyon head 
[ Wallin, 1968]. The author noted some minor changes in the depths of the 
southern end of the subarea over the several dives that were accomplish- 
ed. Significant changes in depth (greater than 25> ft) have been noted 
by Cdr. Don Ferrin, USN (Ret.) [personal communication J several times 
over the past 12 years. The subarea is almost entirely covered with 
coarse sand. Much of this appears to have been derived from outcrops 
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in the near area. 



C. LITTORAL DRIFT 



Wave refraction diagrams constructed indicated that the littoral 
drift in the study area is south. Additionally, the drift along San 
Jose Creek Beach is to the north and a node exists just north of the 
head of the Carmel Canyon. The area of this node was described in a 
previous section (Subarea £)• 

Stream mouths generally tend to migrate in the direction of littoral 
drift [u.S. Army Coastal Engineering Research Center, 1 966 ] . Evidence of 
the southerly drift in the study area was given by a southerly migration 
of the mouth of the Carmel River after opening. The rate of this drift 
was observed to be dependent upon the wave activity and angle of inci- 
dence of the waves upon the beach. Further southward migration of the 
river mouth is prevented by rock outcrops. Evidence of the northerly 
drift along San Jose Creek Beach is given by the long-term northward 
migration of the mouth of San Jose Creek (Fig. 3). 

D . CURRENT STUDIES 

Current data obtained from the Savonius rotor meter was judged to be 
invalid. Moderate to heavy swell persisted during the entire week that 
the meter was in place. The bottom wave surge induced by this swell pro- 
duced an excessive amount of scatter and an unreadable record resulted. 

The ducted meter produced better results than the Savonius rotor 
model. Currents at the chosen station for the period that the meter was 
on the bottom were less than 0.2 knots and this, because of its oscilla- 
tory nature, was judged to be a component of wave surge. Currents in 
water less than 20 ft depth could not be obtained because of the hazards 
to divers and equipment in the surf zone. 
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Even minimal currents can be detected by a diver who is required to 
swim against them. During the entire study with the more than 20 dives 
that were completed, no currents were noted except when entering and 
leaving the surf zone. This observation reinforces the results obtained 
using the ducted meter. 

E. SEDIMENT SIZE ANALYSIS 

The tabulated results as obtained from the computer program are 
presented in Table I. 

1 . Mean Grain Size 

The mean grain size of the samples taken ranged from 2 . 71 0 to 
-1 .580. Six of the samples were classified as fine sand, two as medium 
sand, three as coarse sand and seven as very coarse sand. 

Sand-silt-clay-gravel relationships were computed and plotted 
on a sand- silt -gravel diagram similar to the tetrahedron scheme of 
Krumbein and Sloss [ 1 963 J (Fig.'l5)* "Gravel" refers to particles with 
mean grain sizes coarser than -1 .00 and finer than -6.00. 

In general, mean grain size was larger nearshore and in the 
southern end of the study area (Fig. 16). Data from Carter [l97l] were 
used to assist in the locating of the 30 contour. The seaward bend of 
the mean grain size contours in Subarea 5 probably indicates than sedi- 
ments are being transported seaward in this area. 

2 . Standard Deviation 

The values of standard deviation ranged from 0.^10 to 1 .350* 

One sample was classed as poorly sorted. Two samples were well sorted, 
and 15 samples were moderately sorted. No trends were noted in this 
statistical measure. A mineral analysis of the samples might provide 
an insight into this parameter. 
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SEDIMENT SAMPLE SIZE STATISTICS 
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Gravel 
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Figure 1£. Tertiary Diagram: Sand-Silt-Gravel Relationships. 

(After Krumbein and Sloss, 1963) 




Figure 16. Sediment Mean Grain Size Distribution. (30 contour based 
partially on data from Carter, 1971.) 
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3 . Skewness 

Skewness values ranged from -0.U3 to O. 63 . Three samples were 
classified as very negatively skewed, six as negatively skewed, five as 
nearly symmetrical (one of these was sample Number 18 that was taken in 
the bed of the Carmel River), two were positively skewed and two very 
positively skewed. Sample Number 11, taken across a ripple crest, was 
very negatively skewed (-0.37). Sample Number 12, taken in the trough 
of the ripple from which sample Number 11 was taken, was positively skewed. 
This can be readily understood since positive skewness implies skewed 
toward the fine (more coarse material than fine) and negative skewness 
implies toward the coarse. The wave surge turbulence is greater at the 
crests, thus producing negative skewness. 

U« Kurtosis 

Kurtosis values ranged from 0.91 to 3*21 . One sample was classi- 
fied as extremely leptokurtic, four as very leptokurtic, three as lepto- 
kurtic and eleven as mesokurtic . The samples taken in the southern half 
of the study area and either in shallow water or on the beach face tended 
to be of a leptokurtic nature . 
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IV. SUMMARY AND CONCLUSIONS 



Diver observations in the area showed that the sand transport in the 
study area is generally southward. This was both observed in the surf 
zone and surmised from the evidence of previous sand movement in the 
southern end of Subarea k and the northern end of Subarea £. 

The bathymetric survey provided an indication of sand deprivation in 
Subarea J 4 and a slight seaward bending of depth contours in Subarea 
Comparisons of the author 1 s bathymetric chart with previous charts pro- 
vided no real evidence of slumping or sliding. Previous studies have 
given some evidence for this [Shepard and Emery } 1 9Ul ] • 

Littoral drift was determined by wave refraction diagrams to be south 
in the study area and north along San Jose Creek. A nodal point for sand 
transport exists just north of the head of the Carmel Canyon (Fig. 17). 

Sediment size analysis reinforced the existence of a nodal point for 
sediment transport. 

The two major sources of sedimentary material for the study area ap- 
pear to be the Carmel River and local erosion of granodiorite outcrops. 
The Carmel River provides sediments only during the winter when precipi- 
tation is sufficient to create a flow. Local erosion of outcrops is 
continuous but larger during the winter when winter storms produce in- 
creased wave activity. 

Materials provided by the Carmel River are moved south with the 
littoral drift , mixed with locally derived sediments and deposited along 
the north wall of the Carmel Canyon. Sediments provided by local erosion 
of outcrops south of San Jose Creek Beach are moved north with the 
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littoral drift, mixed with the small amount of sediments provided by 
San Jose Creek, and injected into the Carmel Canyon. 

Winter storms probably provide the necessary energy to induce gravity 
sliding or slumping of deposited material. This material is dumped into 
the head of the Carmel Submarine Canyon for eventual deposition in the 
deep ocean. 
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V. SUGGESTIONS FOR FURTHER STUDIES 



The following studies are presently being conducted within the bay: 

1 . sediment migration within Whaler's Cove (L. Leopold, San Jose 
State College , in progress); 

2. gravity survey of Carmel Bay (A. Souto, NPS, in progress). 

Future beneficial studies of the marine environment of the bay might 

include : 

1 . currents and tagged sand studies; 

2. seismic refraction measurements; 

3. heavy mineral analysis; 

u. carbon, carbonate, and organic nitrogen analysis of sediments; 

!?. a sand budget for the bay. 
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APPENDIX A: SELECTED DIVER LOGS 



DIVERS: Howell, Cronyn DATE: 16 October 1971 

Subareas 1; and % were examined during this dive. 

The divers entered the bay through the surf via a pocket beach ap- 
proximately 300 yards south of the Carmel sewer outfall. Approximately 
50 yards from the beach is an outcrop of granodiorite that is about 20 
yards wide and ranges from 1 0 to i+0 ft in depth. It is heavily encrusted 
with numerous encrusting organisms. The outcrop is lengthy and unbroken 
in its north-south direction and appears to act as a barrier to the sand 
not allowing beach sand to be transported seaward . Seaward of the outcrop 
the sand is noticeably finer and slopes gently downward until a depth of 
60 ft 3 at which depth there was a marked increase in the slope . 

DIVERS: Howell; Millward DATE: 23 October 1971 

Subarea 3 and the southern end of Subarea 2 were examined during this 
dive . 

The divers entered the bay through the surf via a pocket beach ap- 
proximately 300 yards south of the Carmel sewer outfall. The break in 
the slope was located and the divers followed this slope in a northerly 
direction for about 200 yards. A 5>-ft pole was pushed into the sand at 
regular intervals in order to determine sand thickness. In all cases 
the sand cover was greater than 5> ft- The divers then surfaced and swam 
to end of the temporary trestle that was being used to lay the Carmel 
sewer outfall line . Seaward of the trestle the bottom was similar to 
other sand bottoms in the area with small ripple marks and gentle slopes 
with an increase in the slope at approximately 60 ft depth. 
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DIVERS: Howell, Cronyn 



DATE: 29 October 1971 



Subareas 6 and 7 were examined during this dive. 

The divers entered the bay through the surf in the vicinity of San 
Jose Creek. After submerging just south of the wash rock, a very steep 
slope in the sand was noted. Although no angle measuring device was 
used, the slope appeared to be greater than 2S>° . Minor slides were 
readily started with just a kick of a flipper. The western portion of 
the wash rock is heavy granodiorite boulders with a very steep slope. 

No leveling out was noted up to a depth of 120 ft. 

DIVERS: Howell, Cronyn, Millward DATE: 30 October 1971 

Subarea 6 was examined during this dive. 

The divers entered the bay through the surf in the vicinity of San 
Jose Creek. The bottom from the beach to the wash rock was checked for 

i 

the sand channels that were previously noted by Wallin [ 1968 ]. These 
channels were not found nor were any significant sand deposits noted 
directly east of the wash rock. The dive commenced at the same location 
as the dive of 29 October 1971 . A thorough investigation of the western 
side of the outcrop revealed no possible passages for sand to be trans- 
ported. The sand on both sides of the outcrop contained large amounts 
of detritus (crab shelly broken kelp, etc.) and supported a substantial 
tubeworm population (approximate density 20 per ft ). It was concluded 
that most of the sand introduced into the bay must be lost into the 
canyon north of San Jose Creek. 

DIVERS: Howell, Mellor DATE: 30 November 1971 

The divers attempted to launch the boat in Whaler’s Cove, Point Lobos 
State Reserve, but were prevented because of high swell. During the 
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night and early morning these swells had destroyed approximately one 
half of the trestle that had been erected for construction of the 600-ft 
Carmel sewer outfall. Large planks (approximately 6 inches x 8 inches) 
that had been bolted to piles driven into granite were destroyed. Need- 
less to say, a dive was considered inappropriate. It was noted that the 
high tide for that day was an exceptionally high spring tide. This in 
conjunction with the large swell was the reason for the destruction. 

DIVERS: Howell, Mellor DATE: 9 December 1971 

Subarea 1; was examined during this dive. 

The divers entered the bay via the boat ramp at Whaler's Cove at 
Point Lobos State Reserve in the boat. The purpose of the dive was to 
locate a position for placement of current meters. A position was lo- 
cated in 69 ft of water on a sand bottom with a slope of 17 °, approxi- 
mately 300 yards south of the Carmel sewer outfall pier . The bottom was 
very clean with fine white sand and small wave induced ripples. 
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APPENDIX B: COMPUTER PROGRAMS 



COMPUTER DATA CARD FORMATS 

WAVE REFRACTION PROGRAM 
Card 1 : Scale Card 

Col. 1-10 chart scale (FORTRAN FI O.I 4 . format) 
Example : 

^ 833.33 

scale: 833.33 ft = 1 inch 

Card 2: Period and Crest Interval 

Col . 1 -2 wave period in sec 

8-1 1 . crest interval 
Example : 

(To 2 T 0 

period: 10 sec 

crest interval: 2.0 









SEDIMENT SIZE ANALYSIS PROGRAM 



Cards 1 & 2 : Title cards 

Col. 1 -80 of each card contain alphanumeric information to 
appear at the top of the output. 

Card 3 : Identifier for sample 

Col. 1-9 cruise number 
1 0-1 2 sample number 

13-18 sample type Columns 1-28 can contain any 

19-20 depth legal keypunch character. 

21 -22 month 
23-2U day 
25-28 year 



Col. 29-31* latitude (XX XX. XX) 

33-1*0 longitude (XX XX. XX) (First digit of longitude is 

determined from octant . ) 
37-61 depth from top of core(XXXXX.) 

62-66 length of core (XXXXX.) 

79 octant (see below) 

Example : 

f 72 -HOW-01 10DIVER 3o 3231 9723631 7821 3339 l| 



Cruise Number: 


72 -H0W-01 


Sample Number: 


10 


Sample Type: 


DIVER 


■. Depth: 


30 


Date : 


3/23/1 972 


Latitude : 


36°31 .78' 


Longitude : 


1 21 °33- 39 


Sample Detail Cards 





Col. 1*1 -1*1* phi size (absolute value) Fl*.2 (decimal assumed) 
1*3 sign of phi size (+ or -) 

30-36 fraction weight F7«l* (decimal assumed) 

80 end of data flag 

=8 if last phi size for this sample 
=9 if last phi size for all samples 

Example : 

100 + 

Phi Size: +1.0 

Sample Weight : 3 • 9l*36 

Flag: 8 
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Coding for Octant of Geographic Position 



Longitude East 



Longitude West 



o 



0 



1 00°E 



180°W 



100°W 



o 

0 



6 


£ 


1 


2 


oo 


7 


3 


h 



90°N 

Latitude 

North 



latitude 

South 



90°S 



50 



Sample of Computer Output for Wave Refraction Program 



WAVE 


PERIOD = 


11 SECONDS 




WAVE LENGTH 


= 619.52 


FEET 


N = 1 


.00 




J 


SCALE: 1 IN. 


= 833.3 


FEET 


DEPTH 


D/LO 


L/LO 


WAVE LENGTH 


NL 


FM. 


FT. 






ACTUAL 

FT. 


CHART 

IN. 


CHART 

IN. 


1 


6 


0.010 


0.244 


150.86 


0.18 


0.18 


2 


12 


0.019 


0.341 


211.15 


0. 25 


0.25 


3 


18 


0.029 


0.413 


255.90 


0.31 


0.31 


4 


24 


0.039 


0.472 


292.37 


0.35 


0.35 


5 


30 


0.048 


0.522 


323.38 


0.39 


0.39 


6 


36 


0. 058 


0.566 


350.42 


0.42 


0.42 


7 


42 


0.068 


0.604 


374.36 


0.45 


0.45 


8 


48 


6.077 


0 .639 


395.79 


0.47 


0.47 


9 


54 


0.087 


0.670 


415. 12 


0.50 


0.50 


10 


60 


0 .097 


0.698 


432.65 


0. 52 


0.52 


15 


90 


0.145 


0.807 


500.00 


0.60 


0.60 


20 


120 


0.194 


0.878 


543.78 


0.65 


0.65 


25 


150 


0.242 


0.924 


572.17 


0. 69 


0.69 


30 


180 


0. 291 


0.9 53 


590.13 


0.71 


0.71 


40 


240 


0.387 


0.981 


607.69 


0. 73 


0.73 


50 


300 


0.484 


0.990 


613.61 


0.74 


0.74 


60 


360 


0.581 


0.993 


615.46 


0.74 


0.74 



Sample of Computer Output of the Sediment Size Analysis Program. 



MARINE SEDIMENTS NEAR CARMEL RIVER STATE BEACH, CARMEL, CALIFORNIA 

BUFORD F. HOWELL 



CRUISE 72-HOW-Ol SAMPLE NUMBER 5 DEPTH 20 



SAMPLER TYPE DIVER DATE 
DEPTH FROM TOP OF CORE 



5/25/1972 
0. MM. 



LAT. 36-32. 09N 
LENGTH OF CORE 



LONG. 121-55. 72W 
0. MM. 





PHI 


SAMPLE 


FRACTION 


ACCUMULATED 






SIZE 


WEIGHT 


PERCENT 


PERCENT 








-3.00 


1.9983 


4.33 


4.33 








-2. 50 


0. 4966 


1. 08 


5.41 








-2. 00 


1.3831 


3. 00 


3.41 








-1. 50 


3.0940 


6.71 


15.12 








-1. 00 


10. 5151 


22. 81 


37.93 








-0. 50 


11.03 73 


23. 94 


61. 87 








0.0 


4.7591 


10.32 


72.19 








0. 50 


3.2197 


6.98 


79.18 








1.00 


2.0965 


4.55 


83.72 








1. 50 


2.5737 


5.58 


89.31 








2.00 


2. 0364 


4.42 


93.72 








2.50 


1.7854 


3. 87 


97. 60 








3.00 


0.8432 


1.83 


99.43 








3. 50 


0.2150 


0.47 


99.89 


* 






4.00 


0.0499 


0.11 


100.00 








POST-ANALYTICAL 


WEIGHT IS 


46. 1032 






PHI SIZES 


AT PERCENT LEVEL 


OF 








1 C C ) 


5 


16 


25 


50 (M) 


75 34 


95 


-5.21 


-2.68 


-1.46 


-1.21 


-0. 81 


0.17 1.03 


2.15 4 PT. 


-5.21 


-2.63 


-1.47 


— 1.25 


-0.75 


0.19 1.02 


2.13 LINEAR 


SAND, SILT 


, CLAY, 


RELATIONSHIPS 








GRAVEL 


SAND 


SILT 


CLAY 


TOTAL 


SAND/MUD TYPE 




37.93 


62.07 


0 .0 


0.0 


100.00 


9999.99 SAND 




TRASK VALUES 


4 PT . PHIS 








Q1 


Q 2 


Q3 


SO 


LOG SO 


SKG 




2.305 


1.756 


0.887 


1.613 


0.208 


0.81 




INMAN VALUES 












MEDIAN 


MEAN 


DEV. 


SKEW. 


2ND SKEW. KURT. 




-0. 81 


-0.22 


1 .24 


0.48 


0.44 


0.94 




FOLK AND WARD VAL 


UE 










MEAN 


-0 .42 


VERY 


COARSE SAND 






DEVIATION 


1.35 


POORLY SORTED 








SKEWNESS 


0.35 


VERY 


POSITIVELY 


' SKEWED 






KURTOSI S 


1.44 


LEPTOKURTIC 








TRASK VALUES 


LINEAR PHIS 








Cl 


Q2 


03 


SO 


LOG SO 


SKG 




2.384 


1.679 


0.876 


1. 649 


0.217 


0 .86 




INMAN VALUES 












MEDIAN 


MEAN 


DEV. 


SKEW . 


2ND SKEW. KURT. 




-0.75 


-0. 23 


1.25 


0.42 


0. 38 


0.93 





FOLK AND WARD VALUE 



MEAN 
DEVIATION 
SKEkNESS 
KURTCSI S 



-0.40 

1.35 

0 .31 

1 .37 



VERY COARSE SAND 
POORLY SDRTED 
VERY POSITIVELY SKEWED 
LEPTOKJRTIC 
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COMPUTER PROGRAMS 



C 

C THIS PROGRAM IS DESIGNED TO OUTPUT THE NECESSARY 

C INFORMATION REQUIRED FOR MANUAL CONSTRUCTION OF 

C WAVE REFRACTION DIAGRAMS USING THE METHOD GIVEN IN 

C H. 0. PUBLICATION 234. THE ONLY INPUTS REQUIRED 

C ARE: 

C 

C 1 . CHART SCALE 

C 2. WAVE PERIOD 

C 3. CREST INTERVAL 

C 

C THE RATIO L/LO IS PROVIDED BY A SUBROUTINE THAT 

C USES THE NEWTON-S ITERATION METHOD TO SOLVE THE 

C EQUATION L/L0=TANH(2*PI*D/L ) . 

C 



COMMON/ IN PIT/PERIOD, DEPTH 
P E AD (5 , 50 ) SCALE 
• 50 FORMAT (F10.4) 

DO 1000 1=1,5 
READ (5,100) I PE R , ZN 
100 FORMAT ( I2,5X,F4„2) 

PERI 0D = FLO A'T ( I PER) 

C 

COMMENT: CALCULATE THE DEEP WATER WAVE LENGTH (ZLO). 

C 

ZLC=5. 12*IPER**2 

C 

COMMENT: PRINT A TABLE HEADING. 

C 

WRITE (6, 200)1, IPER,ZLC,ZN,SCALE 
200 FORMAT! • 1* ///42X ' TABLE • ,12 ,//15X' WAVE PERIOD = *12, 

*' SECONDS' , 10X, 'WAVE LENGTH = «F7.2,' FEET'//15X'N = • 
* F4 • 2 , 2 6 X , ' SCALE: 1 IN. = • , F6 . 1 , • FEET'/) 

WRITE(6, 205) 

205 FORMAT! 0 '/, 18 X ,' DEPTH • 8X ,' D/L 0 ', 6X, ' L /L 0 ', 9X ,' WAVE • , 
*' LENGTH' ,7X, ' NL* ,// ,15X, 0 FM. ', 3X , 0 FT .', 27X ,' ACTUAL 0 
CHART CHART • /54X, ' FT. IN. IN.') 

K=0 

C 

COMMENT: PROVIDE DEPTHS. 

C 

DO 800 J =1 , 2 00 
210 IFIK.LE.9) GO TO 3CC 
220 IFIK.LE.25) GO TO 320 
230 IFtK.LE.90) GO TO 330 
240 IF(K.LE.IOOO) GO TO 340 
300 K=K+1 

GO TO 4CC 
320 K=K+5 

GO TO 400 
330 K = K + 1 0 

GC TO 400 
340 K=K+25 
400 I FT = K^6 

RAT= I FT/ ZLC 

C 

COMMENT : I F D/LO IS GREATER THAN .6, GO TO THE NEXT TABLE. 

C 

IF (RAT.GT..6) GC TC 1000 
C EPT H= FLOAT ( I FT ) 

C 

COMMEN T : C ALL SUBROUTINE TC CALCULATE L. 

C 

CALL WAVE (WAVED 

C 

COMMENT: CALCULATE L/LO. 

C 

RATIO=WAVEL/ZLG 
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c 

CCMMENT : CALCULATE DATA FOR CONSTRUCTION OF DIAGRAM. 

C 

CHART=WAVEL/SCALE 

ANL=ZN*CHART 

C 

CCMMENT: PRINT RESULTS. 

C 

600 hRITEI 6,70 0 )K, I FT, RAT t RATI C, WAV EL, CHART, ANL 
700 FORMAT ( » 0 • , 14X , I 3 , 3X , I 4 , 6X, F 5 . 3 , 5X , F 5.3 , 6 X ,F7. 2 ,3X , 
*F4.2 ,5X ,F4.2 ) 

800 CONTINUE 

WRITE (6,900) 

900 FORMAT ( ' ' , / , 130 ( 1H-) ) 

1000 CONTINUE 
STOP 
END 



SUBROUTINE WAVE(WAVEL) 

CCN’MCN/I N PUT/ PERIOD, CEPTH 
EXTERNAL FCT 
DDPTH=DEPTH 

CALL ZERO( WAVE L ,FV,DERF ,FCT,DDPTH ,.005,200, IE R) 

RETURN 

END 



SUBROUTINE ZERO(X,F,DERF,FCT,XST,EPS,IEND,IER) 

C 

C PURPOSE 

COMMENT: SOLVE THE GENERAL NCN-L IN EAR EQUATIONS CF THE FORM 
C F ( X ) = 0 BY MEANS OF THE NEWTCN-S ITERATION METHOD. 

C 

I ER=0 
X = XST 
TQL=X 

CALL FCT (TOL , F , DERF ) 

TCLF=100.*EPS 
DO 6 I = 1 , I END 
IF(F.EO.O.) GO TO 7 
IF (DERF. EQ.O.) GO TC 8 
2 CX=F/DERF 
X=X-DX 
TC L = X 

CALL FCT ( TOL , F , DERF ) 

TC L= E PS 
A= ABS (X) 

IFIA.GT.1. ) TOL=TOL*A 

IF ( ABS( DX) .GT.TOL ) GO TC 6 

IF (ABS (F) . LE.T OLF) GO TO 7 

6 CONTINUE 
I ER = 1 

7 RETURN 

8 I ER= 2 
RETURN 
END 



SUBROUTINE FCT ( X , F , CERF ) 

COMMON/I NPLT/PERIOD .DEPTH 
DATA GG/32./,TW0PI/6. 28318/ 

ARG = DEPTH*TWCP I/X 
F= GG--'PERI0D**2/TWCPI-X/TANH(ARG) 

CERF=( -T ANH( ARG J-ARG/COSH ( ARC) *^2) /TANH( ARG) **2 

RETURN 

END 



C** #### # 5|C5*:^C3{C5?:5iC # # #### * # $ $ & & £ $ it 

c 

C SEDIMENT SIZE ANALYSIS MAIN PRGGRAM 

C 

## # :fc ^ ^ # 3$c 5jc ajc 3fc ?{c s}s 2$c s|e :fc £ # ;{c sje sjc # ak # # # £ ^ 

c 



COMMON /BLK1/PHIS( 8,2) , N I NM , NF A Li , N SSC , SUMN L ,NTR SK 
REAL NEG/'- •/ 

DATA ANORTH/’N '/, SOUTH/’S '/,EAST/'E '/, 
*WEST/'W •/ 

DIMENSION FRWTUOO) ,PRCT(100 ) , CRMC (20 ) , STMC ( 20 ) , 

*EXMC (20) ,CRMB( 20) ,STMB(2 0) ,EXMB(20) ,ICR(20),JCR(20), 
*TATLE(20 ) ,T ITL E { 20 ) 

REAL*8 CRUZ ,CRUZ R ,CRMB , CR MC , S MPLR 

DIMENSION T(IOO) , PH I ( 1 CO) , AC PC ( 1 00 ) 

DATA KKK/O / , MAT / 0/ , MBT/O/ ,MED/0/ , KSM/O/ 

NI NOT = 1 

902 N2= 1 
NA=1 
NZ = 1 
NL = 1 
KB = 1 
M= 0 

PAWT=0. 

DO 903 K=I *100 
T ( K ) =4. 09 
ACPC { K ) =999 .99 
FR WT { K ) =0. 0 

903 PH I ( K ) = 99. 99 
MA=0 

MB = 0 

MBT=0 

MC=0 

MED=0 

SUMWT=0 .0 

SUMPC=0 .0 

AC PC { 1) =0. 0 

R E AD { 5 , 20 ) TITLE 

20 FORMAT ( 20A4 ) 

R E AD { 5 , 2 1 ) TATLE 

21 FORMAT (20A4) 

1 f\K = l 

2 READ { 5, 26)CPUZR ,MCR, STATR ,SMPLR ,EXID , MO , D A ,YR , LAT A, 
*CEGLT , LNGA,DEGLN,PHIR, S 1 GN , B FWT , DFPTHR , CP LN , IQUD,NE 

26 FORMAT (A7»A2»A3,A6»3A2, A4 ,I2,F4.2» 12 , 2F4. 2 , A1 , 4X, F7.4, 
*2 F 5 . 0 , 12X, 211) 

T IMES = 1 . 

I F ( SIGN. EQ. NEG) TIMES=-1. 

PHIR=(TIMES*PHIR) 

3 IF (PHIR.LT.3.9) GO TO 29 
IF(NZ.EQ.l) GO TO 48 
IF(PHIR.LE.4.1) GO TO 46 
NZ=2 

GO TO 29 
46 NZ=1 

GO TO 44 

48 IF(NL . EQ.2 ) GO TO 44 
N L =2 
NZ = 2 

29 FRWTR=B FWT 

I F ( FRWT R .GE.O . ) GO TO 52 
I F ( (FRWTR+.C1) .GE.O.) GO TO 25 
WRITE(6,51) FRWTR 
25 F RWTR=0. 0 
52 PAWT=PAWT+FF.WTR 
SU MW T = P A WT 
GO TO 49 

44 I FINK. EQ.2) GO TO 59 
NK = 2 

FWT1=BFWT*50. 

GO TO 2 

59 FWT 2= B FWT *50 . 
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FRWTR=FWT1-FWT2 
I F ( FRWTR .GE . 0. ) GO TO 53 
I F ( ( FRWTR+O .01 ) .GE.O. ) GO TO 55 
60 WR I TE ( 6 , 51 ) FRWTR 

51 FORMAT ( 25X , 'WEIGHING ERROR. FRACTION WEIGHT WAS ', 
*F10.4,' BEFORE BEING SET TO ZERO') 

£5 FR WTR = 0 . 0 

53 FWT1=F WT2 
NI NOT =2 

IFCNE.LT. 8) GO TG 52 

54 PAWT = PAWT + FRWT R+FWT2 
SUMWT = P AWT 

K = K+ 1 

M=M+1 

KK=KK+ 1 

PHI (K)=PHIR 

FRWT ( K )=FRWTR 

K = K+1 

M=M+ 1 

KK = KK + 1 

FRWT(K) =FWT2 

PHI (K ) =12. 

IFCPHIR.LE.il.) PHI (K)=PHIR+1. 

GO TO 100 

49 I F ( NA . EQ .2 ). GO TO 50 
NA=2 

C DETERMINE OCTANT FROM FIRST CARD OF SAMPLE 

DK = WE ST 
I DH=LNGA 

I F ( I QUD-2 ) 8,9,11 

8 IDH=I DH+100 

9 DG= ANORT H 
GO TO 23 

11 I F ( I QUD— 4) 12,13,15 

12 I DH= I DH+ 10 0 

13 DG=SQUTH 
GO TO 23 

15 D K = E A S T 

I F ( I QUD-6) 17,16,19 
17 I DH= I DH+ 100 

16 CG=ANORTh 
GO TO 23 

19 IFC IQUD.GE.8) IDH=IDH+100 
DG = SOUT H 
23 CRUZ-CRUZR 
NCR=MCR 
ST AT = ST ATR 
EXC=EX I D 
DP TH=DE P TH R 
WRITE (6,799) TITLE 

799 FORMAT ( • 1 ' , 2 6X , 2 0 A4 ) 

WR I T c ( 6 , 800) TATLE 

800 FORMAT CO* ,26X,20A4//) 

WRITE (6,801) CRUZ, NCR, STAT, EX C , $ MPLR , MO , D A , YR, L AT A , 
*DEGLT,DG, I DH , D EGLN , DK , DE PTHR , CR LM 

801 FORMAT {/36X, 'CRU ISE ',A7,A2, • SAMPLE NUMBER ',A3, 

*• DEPTH SA2//25X,' SAMPLER TYPE *,A6,' DATE ',A2,'/', 
*A2 , ' / ' , A4, ' LAT. • , I 2 , •- ' , F 5 . 2 , A 1 , • LONG. ',13. 

,F5.2,A1/26X, 'DEPTH FROM TOP OF CORE ',F7.0, 

*' MM. LENGTH OF CORE ',F7.0,' MM.') 

NB = 2 
KK=0 
K = 0 

IF ( FRWTR .NE.O. ) N2=3 
PH I A=-12 .0 
GO TO 1 

C SAVE PHI 8 FRACTION WT . FROM EACH DETAIL CARD 

50 KK=KK+ 1 
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K =K + 1 

PHI ( K) = PH I R 

FRWT ( K ) =FRWTR 

IF (NE.G E.8 ) GO TO 100 

C SET FLAG (N2=2) IF NEW PHI LESS THAN LAST PHI 

IF(PHIA.GE.PHI R) N2=2 
P H I A=P H I R 
M = M + 1 

IF (NINOT-1 ) 1,1,2 
ICO M = K 

C DETERMINE FRACTION PERCENTS AND ACCUMULATED * 

DO 67 K= 1, M 

PRCT (K ) = ( F RWT ( K ) /PAWT ) *100 . 

SUMPC=SUMPC+PRCT(K) 

ACPC ( K ) =SUMPC 
67 CONTINUE 

CALL TVAL(T,ACPC,M) 

IFISUMPC.LT .99 .94.0R.SUMPC.GT .100.06) GO TO 1500 
I F ( M2-2 ) 103,1666,1990 

103 WRITE { 6 , 802) ( PH I ( J ) , FRWT ( J ) , PRCT ( J ) , ACP C ( J ) , J=1 , KK ) 
802 FORMAT (• O' ,40X ,' PHI SAMPLE FRACTION ACCUMULATED’, 
*/41X,'SIZE WEIGHT PERCENT PERCENT' / (40X , F6 .2 , 

*F8.4,F8.2,F12. 2) ) 

K JK=KK 

WRITE (6,79) PAWT 

79 FORMAT ( /40X, ' POST-ANALYTICAL WEIGHT IS’,F9.4) 
SUMNL=ACPC(KK) 

N S SC=2 
NTRSK= 1 
N I N M = 1 
NF AW= 1 
JJ = 0 

DC 77 I=1,KK 

I F ( FR WT ( I) . EQ. 0. ) GC TO 77 
J J=J J+l 
T(JJ)=T (I ) 

ACPC ( J J ) =ACPC ( I ) 

PHI( JJ )=PHI ( I ) 

77 CONTINUE 
K K = J J 
K JK= J J 

IF (SUMNL.GE.72.0) GO TO 105 
N SSC= 1 

W R IT E ( 6 , 809 ) PHI(KK) 

839 FORMAT (' 0’ ,29X ,' DID NOT INTERPOLATE ANY PHI SIZES', 

*' BECAUSE' /29X, ' ACCUMULATED PERCENT AT ',F5.2, 

*' DID NOT EXCEED'/29X, '72 PERCENT') 

GO TO 163 

105 IFI4.LT. KK) GO TO 111 
WRITE(6,815) KK 

815 FORMAT ( ' ',28X,' ONLY ',13,' DETAIL CARDS SO ONLY', 

*' SAND. '/30X, ' SILT , CLAY RELATIONSHIPS CALCULATED.') 

GC TO 901 

C COMPUTE T-VALUES 

111 CALL I NTP P ( ACP C , PH I ,T , KK ) 

160 CALL SNSTCLtPHI , ACPC , SUMPC , KJK ) 

CALL CTIFW 

C INCREMENT COUNTERS AND GO TO NEXT SAMPLE 

931 KFK=KKK+KK+1 
K SM = K SM+ 1 
1 = 1 + 1 
N I NOT = 2 

IFINE.NE.9) GO TO 9C2 
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C PRINT RESULTS AND MESSAGES 

950 WRITE(6,860) KSM.KKK 

860 FORMAT {' 1' ,30X ,' THIS BATCH OF CARDS CONTAINED DATA', 

*• FROM ',14,' SAMPLE S ' /30X t ' FOR A TOTAL OF ',14, 

CARDS.' ) 

IF(MST.GT.O) WRITE(6 t 863) ( C RM B { L ) t I CR ( L ) t STMB ( L ) , 
*EXMB(L ) ,L=1,MBT ) 

863 FORMAT {• O' t 29X CARDS OUT OF ORDER ON THE FOLLOWING' , 
1* SAMPLES'/31Xt 'CRUI SE' , 5X , ' SAMPLE ', 5X , ' DEPTH' // t 

2 (33X, A7, A2,A3, A2 ) ) 

IF(MC.GT.O) WR I T E ( 6 , 864 ) (CRMC(L) , JCR ( L ) , STMC ( L ) , 
lEXMC(L) , L= 1 , MC ) 

864 FORMAT {• O' , 29X , 'NO ZERO PERCENT CARDS ON THE', 

1' FOLLOWING STATIONS' /31X, 'CRUISE, 5X , 'SAMPLE ' , 5X , 

2 'DEPTH' //( 33X,A7,A2,A3,A2) ) 

MTT=MBT +MC+MED 
IF(MTT.GT.O) GO TO 959 
WRITE(6,861) 

861 FORMAT (' O' ,30X ,' CONGRATULATIONS NO ERRORS WERE', 

1' FOUND IN THIS BATCH OF CARDS') 

GC TO 960 

959 WRITE (6,865) MTT 

865 FORMAT ( 'O' ,29X, ' SORRY OLD CHAP, BUT YOU MADE ',13, 

1* ERRORS ON THE DAT A • /30X , FOR THIS RUN. NEXT TIME BE ' 
2' MORE CAREFUL') 

960 RETURN 

1666 WRITE (6 ,8666 ) 

8666 FORMAT! • O' ,29X ,' CARDS CUT CF ORDER. CHECK VALUES', 

1 'BELOW') 

MBT=MBT+1 
CRMB ( MBT ) =CRUZ 
I CR ( MBT ) = NC R 
STMB ( MBT ) = ST AT 
EXMB (MBT ) = E XC 
C-0 TO 1501 
1990 WR I TE ( 6,8990) 

8990 FORMAT ( 'O' , 28X , ' NO ZERO PERCENT C ARD . • /3 OX , • CHEC K • , 

1 • VALUES BELOW. • ) 

MC =MC+ 1 
CRMC ( MC ) =CRUZ 
JCR( MC) =NCR 
STMC ( MC ) =STAT 
EXMC( MC ) = EXC 
GC TO 1501 

1 5C0 WRITE(6, 830 PAWT 

830 FORMAT (' O' , 32X , 'SUM CF FRACTION WEIGHTS DID NOT ', 

1 ' EQUAL POST ANALYTICAL WEIGHT '/32X, 'WHICH WAS ', 

2F8 .3 , ' . CHECK THE VALUES BELOW FOR ERRORS.') 

1501 WR ITE (6 , 83 1 ) 

831 FORMAT ( ' O' ,40X , ' PHI FRACTION FRACTION 4CCUM. T-', 
1/41X,' SIZE WEIGHT PERCENT PRCT VALUE') 

WRITE (6,833) (PHI(J), FRWT ( J ) ,PRCT ( J ) , AC PC ( J ) , T ( J ) , 

1 J=1,KK) 

833 FORMAT CO* , 40X , F 5 . 2 , F9 . 3 , F9 . 2 , F 8 . 2 ,F 7 . 3 ) 

WRI TE ( 6 , 8333 ) SUMWT 

8323 FORMAT ( 'O', 29X , ' SUM FRACTION WEIGHTS = ',F8.3,' GRAMS' 
1 ) 

MED=MED+ 1 
GO TO 901 
999 RETURN 
END 
C 
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INTERPOLATION SUBPROGRAM 



SUBROUTINE I NT R P ( AC PC , PH I , T , KK ) 

DIMENSION ACPC< 1) »PHI (1) ,T(1) 

COMMON /BLK1/PHIS(8,2) ,N I NM , NFA W , N SSC , SUMNL ,NTRS K 
REAL PC (8 ) /I ., 5. ,16. ,25 .,5 3 .,7 5., 34., 95./, 

1 XTT ( 3)/-2. 3 25, -1.645, -.99 5, -.674, 0.0,. 674.. 995, 1.645/ 
PHIS (6, 1 ) = 99 .99 
PHIS (6,2)=99.99 
PHIS (7, 1 ) = 9 9 . 9 9 
PHIS(7,2) = 99.99 
PHIS(8,1)=99.99 
PHIS (8,2)=99.99 
N I NT=0 

1101 IF(NINT.EQ.O) GO TO 1102 
PHIS (NINT, 1)=YYA 

PHIS (NINT,2)=YYL 

1102 N I NT=N I NT+ 1 
XPC=PC( NINT ) 

XT = XTT (NINT.) 

DO 151 L=1,RK 

IF (ACPC( L )-XPC ) 151,152,153 

151 CONTINUE 

152 YYA=PH I ( L ) 

YY L = PH 1 ( L) 

GO TO 210 

153 1 F ( X PC-7 5.0) 154, 155, 157 

155 IF(SUMML.GE.75.0) GO TO 156 
N 1 NT=9 

GO TO 1197 

157 IF (XPC.GT.84.0 ) GO TO 159 
I F( SUMNL.LT. 84. ) GO TO 160 

156 IF(L-KK) 195,197,197 

159 I F ( SUMNL-95. 0) 169,156,156 



168 


I F ( SUMNL.LT. 81 
N I NT=1 0 
GO TO 1197 


.0) 


GO 


TO 


1475 


169 


I F ( SUMNL .LT .92 
M NT = 1 1 


.0) 


GO 


TO 


1476 


1197 


XT = XTT( NINT-3) 
GO TO 197 










154 


IF(L.LE.2) GO 


TO 


196 






195 


LA=2 

GO TO 199 










196 


LA= 1 

GO TO 199 










197 


L A = 3 










199 


LS=L-L A 











X=XT 

X 1 =T( LS) 

X2=T( LS+1) 

X3 = T ( LS + 2 ) 

X4 = T( L S + 3 ) 

Y 1 = PH I ( LS ) 

Y2 = PH I ( LS+1 ) 

Y 3 = P H I (LS+2) 

Y4 = PH I ( LS +3 ) 

C AITKENS FOUR POINT INTERPOLATION 

P12 = ( ( Y 1* ( X2-X ) ) - (Y2* ( Xl-X ) ) ) / ( X2— X 1 ) 

P 13= ( ( Y 1* ( X3-X) )-( Y3*( Xl-X) ) ) /( X3-X1 ) 

P 14= ( ( Y 1* ( X 4— X ) )— ( Y4* ( Xl-X) ) ) /(X4-X1) 

P 1 2 3 = ( ( PI 2* ( X3-X ) )-(P13* (X2-X ) ) )/ (X3-X2 ) 

P 124= ( ( P 12* ( X4-X ) )-(P14* (X2-X) ) ) /( X4-X2) 
YY A= ( ( P123* (X4-X) )- ( P124*( X3-X ) ) )/(X4-X3) 
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C LINEAR INTERPOLATION 

XI 1=T ( L — 1 ) 

X22=T ( L ) 

Y11=PHI (L-l) 

Y22=PHI ( L) 

YYL=( X-X11)*( Y22-Y11) /( X22-X11 )+Yll 
IF(L.LE.2) YYA= YYL 
210 IF (N INT. LT. 8) GO TO 1101 
IN I=N INT— 7 
PHIS (NINT, 1 }=YY A 
PHIS(NINT,2)=YYL 

WRITE (6, 803) ( ( PHI S ( 1 1 , K ) , 1 1 =1 , 8 ) , K= 1 , 2 ) 

803 FORMAT (' O' ,29X ,' PHI SIZES AT PERCENT LEVEL OF ' 

1 /32X,'l (C) 5 16' ,7X,'25' ,7X,'50 (M) 75' 

2,7X, '84', 7X, '95'/27X,8F9.2, ' 4 PT. 

3' LINEAR') 

RETURN 

1275 N M3=N I NT-3 

PHIS (NM3,1 )=YYA 
PHIS (NM3,2)=YYL 
I PH=NM3 
I PE R=P C ( NM3 ) 

IF(NINT.EQ.ll) N FA W = 2 

1300 WRITE (6,1301) I PE R , ( ( PH I S ( 1 1 , K ) , 1 1= 1 , 8 ),K=1,2) 

1301 FORMAT ( ' O' ,29X , ' PHI SIZES AT PERCENT LEVEL OF (',12, 
1' LEVEL EXTRAPOLATED) '/32X, 1(C) ' ,6X,' 5' ,7X,' 16' ,7X, 

2' 25' , 7 X , ' 5 0 ( M ) 75 ' ,7X, • 84' /30X , F5 .2 , 7F9 . 2, • 4 PT.'/, 

330X,F5.2,7F9.2, ' LINEAR') 

RETURN 

1475 I P E R = 8 4 
N I N M = 2 

GO TO 1400 

1476 NF AW=2 
IP EP.= 95 

GO TO 1400 

1400 WRITE (6,1401) I PE R , ( ( PHI S ( I I , K) , 1 1 =1 , 8 ),K=1,2) 

1401 FORMAT (' O' , 29X ,' PHI SIZES AT PERCENT LEVEL OF (',12, 

1' LEVEL EXTRAPOLATED) '/32X,1 (C )' ,6X, '5' ,7X, • 16', 7X, 
2'25',7X,'50(M) 7 5 ' , 7 X , ' 84 « /3 OX , F5 . 2 , 7F9 . 2 , • 4 PT.'/, 

330X, F5.2,7F9 .2, ' LINEAR') 

RETURN 

END 



C 

C nC JU «X» %JU iX- X X X X X X X X X# X X X X X X X X X X X x x X X X X X X ■'X x x X x x x x X X X x X* x X 
-T* n * - % -v* ^ # r -f ^ t ^ ^ ^ t t ^ ^ ^ ^ t t t ^ t t ^ t ^ -p v t ^ t t ^ 

c 

C SUBPROGRAM WHICH IS THE EQUIVALENT OF PLOTTING 

C HAND, THE GRAIN SIZE AGAINST THE ACCUMULATED 

C PERCENTAGE. 



n* **5^ 



, BY 



X X* w x* X x % 

-r -r nr -V - 



SUBROUTINE TVA L ( T, ACPC ,M) 

CIMENS ICN T ( 1) , ACPC( 1) 

REAL TBL PC ( 87) / 0.0, 01. 20, 03. 59, 04 .78, 05.96, 07. 14, 00.32 

* ,09.48, 

1 10. 64, 11. 79, 12. 93, 14. 06, 15. 17, 16. 28, 17. 36, 18.44, 19. 50 

* 20.54, 

2 21. 57, 22. 57, 23. 57, 24. 54, 25. 49, 26. 42, 27. 34, 28. 23, 29. 10 

* ,29.95, 

3 3 0. 78 ,31. 59 ,32. 3 8,33. 15 ,33. 89 ,34.61 ,35 .31 ,35 .99 ,36 .65 
*, 37 . 29 , 

4 37 .90 ,38.49,39 .07 , 39. 62, 40. 15, 40. 66, 41. 15,41.62, 42. 07 

* ,42.51, 

5 42.92, 43.32,43.70,44. 06,44. 41,44. 74,45. 05,45.35,45.64 

- ,45.91, 

6 46. 16,46.41,46.64,46. 86,47. 06,47.26,47.44,47.61 ,47.78 

* ,47 .93, 

7 4 8. 08 ,48.21,48.34 ,48.46,48. 57,48.68 ,48.78,48. 87,48. 96 

* ,49.04, 

8 4 9.1 1 ,49 .20,49 .31 ,49.4 0, 49. 60,49. 70, 49. 8 0, 4 9.90, 50. 0/ 
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REAL TBLT< 87)/ 0.00, 0.03, 0.09, 0.12 , 0.15,0.18, 0.21, 0.2*, 

* 0.27,0.30, 

1 0.33 ,0.36,0.39, 0.42 ,0.45, 0.48, 0.51, G.54,0.57, 0.60, 

* 0.63,0.66, 

2 0.69,0.72, 0. 75, 0. 78, 0. 81 , 0. 84,0. 87,0.90,0.93 ,0.96, 

* 0.99,1.02, 

3 1. 05, l.C 8, 1. 11, 1.14, 1.17, 1.20, 1.23, 1.26, 1.29, 1.32, 

* 1.35,1.38, 

4 1.41 ,1.44, 1.47 ,1.50,1. 53 ,1.56,1 .59,1 .62,1 .65, 1 .68, 

* 1.71,1.74, 

5 1.77,1.80, 1.83,1.86,1.89,1.92, 1.95, 1.98,2.01,2.04, 

* 2.07,2.10, 

6 2. 13, 2. 16, 2. 19, 2. 22, 2. 25, 2. 28, 2. 31, 2. 34 ,2. 37, 2. 41, 

* 2.46,2,51, 

7 2.65,2.75,2.88,3.08,4.09/ 

C CALCULATE T-VALUE 

CO 68 K= 1 , M 

IF (ACPC(K) .LT. 100.0) GO TO 62 

61 T { K ) =4 . 09 
C-Q TO 68 

62 DLPC=ACPC( KJ-50. 

IF(DLPC) 64,63,65 

63 T { K ) =0 . 0 
GO TO 68 

64 GMPC=— DL PC 
GC TO 66 

65 GMPC =D LPC 

66 DO 69 L= 1 , 8 7 
IF(TBLPC(L )— GMPC ) 69,70,71 

69 CONTINUE 

70 TCALC=TBLT { L ) 

GC TO 72 

71 TCALC=< GMPC- TB LPC ( L-l ) ) * (TBLT < L ) -TBLT (L- 1 ) )/ 

* TBLPCIL )-TBLPC(L-l) )+TBLT(L-l) 

IFtTCALC.GT .4.09) GO TO 61 

72 IF(DLPC.GE.O.O) GO TO 74 
TCALC=-T CALC 

74 T { K) =T CALC 
68 CONTINUE 
RETURN 
END 
C 

C -x x x x x x x x x x -x> x x x x x x x x X x x X x x X x X X x X x x x x X x x x X - X x x x x x x x x x x x x 

'V' ■*,' «v »■,*■ ■*> X- ■<-,> x- 'E •T’* x 1 ^ 9 nr* O' X ^ X 'i' v X* *r X 1 'r 'T 4' '' “V* 'I s ■X V n' *r ■X "r X *p X X V t T ^ t t t 'f -V •X "V* n* 

C 

C SUBPROGRAM FOR COMPUTING SANC-S ILT-CLAY RELATIONSHIPS 

C 

/"* *sJL» st. sJU vlU vb *JU si# J# si# sV* 4L# si# 4/ si# 4# 4# 4U 4# si# si# si# JU s’# »l# J# J# si# «JU <4# si# 4# si# O# sJL 4# >JL# si# 4^ 4# *X* «J#» 4# O- JL> 4# 4# si# 4# 4# 4# J# 

I •*■/** #|S #-ys #,s #js #|S #p #^s #lp* #JS #|S #p #|S #p ^s #|S #jS rj' #p #p #p #p #p #js #p #p #j» #p *ps As #jS #p #p #p #p «#j% #p Jp #^*» #)S rp #p <+f- #ys #p 



SUBROUTINE SNSTCLI PHI , ACP C , SUMPC , K JK ) 

DIMENSION PHI (1 ) ,T(1) ,ACPC(1 ) 

REAL*8 CLASS! 9) /’ SAND SILT CLAY SANDY 1 
1, 'SILTY CLAYEY -SILT -CLAY •/ 

INTEGER SUBS{33)/1 ,7,7 ,2,7 ,7 ,3, 7, 7, 1,8, 9, 5, 1,7,4, 2,7, 6 
*,1,7, 6, 2, 7, 

14,3,7,5,3,7,7, 7,7/ 

N I =0 

RAT 10 = 9 999. 99 
GRSN=0 .0 
S AND = 0 . 0 
S I L T = 0 . 0 
C L AY=0 .0 
SANDP=0. 0 

IF(PHHKJK) .LT.-1.0) GO TO 380 
CO 302 KG=1,KJK 

IF (ABS ( PHI { KG) +1 .0 ) .LE . .0001 ) GO TO 303 
IF{PH1 < KG) +1.0 ) 302,303,307 
3 02 CONTINUE 
303 GR SN=AC PC ( KG ) 

307 IF (PHI (KJK) . LT .4.0) GO TO 312 
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CC 309 KS=1,KJK 

I F ( AB S ( PH I ( KS ) — 4 • 0 ) .LE..0001) GO TO 310 
IFtPHI (KS)-4.0) 309,310,313 

309 CC NT I MU E 

310 SAND=AC PC ( KS ) 

SANOP=SAND-GRSN 
FMUD=SUMPC-S AND 

I F ( FMUD .NE. 0. ) RATI 0=SAND/FMUD 
GO TO 313 

312 SAND=SUMPC-GRSN 
S ANDP= SAND 

GO TO 380 

313 IF(ABS(PHI(KJK)-8.0).LE«.0001) GO TO 314 
IF(PHKKJK) .LT.8.0) GO TO 317 

314 CO 315 KSL = 1 ,KJK 

I F ( ABS ( PHI ( KSL )-8. 0) .LE..0001) GO TO 316 
IF(PHI(KSL).GE.8.0) GO TO 316 

315 CONTINUE 

316 SILT=ACPC( KSL) -SAND 
CLAY=SUMPC-SAND-SILT 
GO TO 320 

217 SILT=SUMPC-SAND 
320 N 1 = 1 

IF (SAND. GE. 75.0) GO TO 380 
N I =2 

I F ( S ILT .GE .75 .0 ) GO TO 380 
NI =3 

IF(CLAY.GE.75.C) GO TO 380 
N I =4 

I F ( S AND o GE . 2 0. 0. CR . SI LT. GE . 2 0 . 0 . OR . CL AY . G E .20 . 0 ) 

* GO TO 380 
N4 = l 

IF(CLAY/SILT.LT.l.) N4=2 
IF(SAND/SILT.LT.l. ) GO TO ( 34 1, 334 ) , N4 
GO TO (334,336) , N4 

334 IF(CLAY/SAND. LT. 1. ) GO TO (338, 337), N4 
GC TO (340, 339), N4 

336 NI=5 

GO TO 380 

337 N I =6 

GC TO 380 

338 N I =7 

GO TO 380 

339 N I =8 

GO TO 380 

340 N I =9 

GO TO 380 

341 N I = 1 0 

380 KPITE (6 ,804) 

804 FORMATt ’ 0* ,28X, • SAND, SILT, CLAY, RELATIONSHIPS' /30X, 
1' GRAVEL SAND SILT CLAY TOTAL 1 

2,' SAND/MUD TYPE') 

NN=3* ( N I- 1 ) + 1 
N1=SUBS ( NN ) 

N2 = SUB S ( NN+1 ) 

N3= SUBS ( NN + 2 ) 

WR ITE (6 ,87 ) GR S N , S AN'DP , S ILT, CL A Y , SUMPC , R AT 10 , CLA S S ( N 1 ) 

* CLA SS ( M2) , 

1 CLASS (N3) 

87 FORMAT ( ' • , 29X , F8 . 2 , 5F9 .2 , 3X , 3 A6 ) 

RETURN 

END 
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OOOOO 



SUBPRCGRAM F CR CCMPUTING TRASK, INMAN, AND FOLK AND 
WARD STATISTICS. 



SUBROUTINE CTIFW 

CCMMCN /BLK1/PHIS (8,2 ) , N INM, NF AW , NSSC , SUMNL , NTRSK 
REAL* 8 ZMEN ( 20) /' VERY COARSE SANDCCARS5 SAND MECIU 

IN SAND FINE SAND VERY FINE SAND COARSE SILT 

2 MEDIUM SILT FINE SILT VERY FINE SILT 

3CLAY •/ 

RE AL*8 DEV 1(18 )/'VFRY WELL SORTED WELL SORTED 

MODERATELY SORTED POORLY SORTED 

VERY POORLY SORTED EXTREMELY POORLY SORTED'/ 

RE AL*8 SKEW (15 )/ 'VERY NEGATIVELY SKEWED NEGATIVELY SK 
1EWED NEARLY SYMMETRICAL POSITIVELY SKEWED 

2 VERY POSITIVELY SKEWED '/ 

RE AL*8 KURT (15 ) /' PLAYTKURT IC MESOKUPTIC 

1 LEP TCKURTI C VERY LE PTCKURT I C 

2 EXTREMELY LEPTOKURTIC '/ 

REAL TLFC(7)/0. 0,0. 35, 0.50,1. ,2 . ,4. ,99.99/ ,TLFS( 6 ) /-3. 

^ -p 

1 I3^9l99/’,TLFK(6 )/0., .9, 1.11, 1 . 5 , 3 . , 99 . 99/ 

REAL HEAD( 6 ) / ' 4 PT. PHIS LINEAR PHIS •/ 

INTEGER ITLFSC(6)/2*1,2,3,4,5/ 

DC 1000 LL=1,2 
IF(NSSC.EQ.l) GC TC 400 

C CALCULATE TRASK VALUES 

PF I 5=PH I S ( 2 , LL ) 

PHI16=PH IS (3 , L L ) 

PHI25=PHIS(4,LL) 

PHI50=PHIS( 5,LL) 

PF 175 = PH IS ( 6 , L L ) 

PHI 84 = PH I S ( 7 , L L) 

PHI95=PHIS(8,LL ) 

Ql=2.** (-PH125 ) 

Q2=2.**(— PHI 50) 

G3=2.**(-PHI75 ) 

SC=SQRT (Q1/Q3) 

FLGS0=AL0G1C( SO) 

SKG=SQRT ( ( Q1*Q3 ) / ( Q2*Q2 ) ) 

L L L =3* ( L L- 1 ) + 1 
LLT=LLL+2 

WRITE (6 ,805) (HEAD(LLK) , L LK = L LL , LLT ) , Q 1 , Q 2 , Q3 , SO, F LGSO, 
* 29 X, ' Q 1 Q2 Q3 SC LOG SO 

*SKG'/30X,F6.3, 4F9.3,F8.2) 

IF (NTRSK. GT.l ) GC TC 1000 

C CALCULATE INMAN VALUES 

IF (NINM.GT. 1) GO TO 505 
F I MD= ( PHI 16 + PH I 84 ) /2.0 
FIDV=(PHI84-PHI16)/2.0 
FISK=( FIMD-PHI 50) /FID V 
I p ( PHI95 .EQ . 99.99) GO TO 504 
F 2 SK= ( ( (PH 1 9 5 + PH 15 )/2. ) -PH I 5 0 ) / F I DV 
FIKU=( ( (PHI95-PHI5) /2. J-FIDV J/FIDV 
WRITE (6, 806 ) P H 1 50 , F IM D , F I DV , F ISK , F2SK, F I KU 
806 FORMAT ( • O' ,29X , ' I NNAN VALUE S ' /32X , ' MED I AN MEAN',5X, 
l’CEV. SKEWED. 2ND SKEW. KURT.'/32X, 

2F5.2,3F9.2 ,2F1 1 .2 ) 

GC TO 600 

504 WRITE(6,92) PH I 50, F IMD » F I D V , F I SK 

1 ' NO SKEWNESS AND KURTOS IS )' /30X ,' MEDIAN MEAN ' 

1, 'DEV. SKEW. ' /30X,F5.2 ,3F9.2) 

GC TO 600 

505 WR I TE ( 6 , 93 ) 
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93 FORMAT { ' O' , 28X, ' INMAN PLUS FOLK AND WARD VALUES NOT • 

1 ' CALCULATED BE C AUS E • /29X , • NEXT TO LAST ACCUMULATED ' 

2PERCENT WAS LESS THAN 84') 

GO TO 1000 
400 WP ITE (6,91) 

91 FORMAT! • O' ,28X, • NCT ABLE TO CALCULATE TRASK, INMAN, • 
*,' OR FCLK AND WARC/29X, 

* • VALUES BECAUSE NEXT TO LAST ACCUMULATED PERCENT DID* 
*,' NOT ' /29X, ' EXCEED 72') 

GO TO 1000 

C CALCULATE FOLK AND WARD VALUES 

600 IF(NFAW.EQ.l) GO TC 601 
WRITE (6,94) 

94 FORMAT (' O' , 29X ,' COULD NOT CALCULATE FOLK AND WARD' 

*' VALUES BECAUSE NEXT' /30X , ' TO LAST ACCUMULATED PERCEN 
*, 'NT DID NOT EXCEED 92. ' ) 

C COMPUTE MEAN AND DETERMINE CATEGORY 

601 FMZ = ( PHI 16+PHI50+PHI84) /3.0 
DC 300 1=1,10 

ZV AL = I— 1 

IF( FMZ.LT.ZVAL ) GO TO 301 

300 CONTINUE 
1=10 

301 11=2*1 I-D+l 
12=11+1 

WRITE ( 6 , 8 C7 ) FMZ,(ZMEN(II) ,11=11,12) 

8C7 FO RM AT ( ' 0 ' , 29X , 'FOLK AND WARD V ALUE ' /3 IX , 'MEAN • , F 12 . 2 , 
1 4X , 2A 8 ) 

C CGMPUTE DEVIATION AND DETERMINE CATEGORY 

FDEV = ( PH 18 A- PH I 16) /4 . + ( PH I 95-PH I 5 ) /6 . 6 
DO 604 L=1 , 7 

IF{TLFD( L) .GE.FDEV) GO TO 605 

604 CONTINUE 

605 IFDTL=L-1 

I F ( IFDTL.EQ.O) I FDTL=1 
1 1=3*( IFDTL-1) +1 
12=11+2 

WRITE (6,131) FDEV, (DEVI (II) ,11=11,12) 

131 FORMAT (3 IX, « DE V I AT I ON • , F 7 . 2 , 4X , 3 A 8 ) 

C COMMUTE SKEWNESS AND DETERMINE CATEGORY 

FSK=(PHI16+PHI84-2 .0*PHI50)/( 2.*( PHI 84-PH I 16) ) + 

1 ((PHI 5+PHI95 )-(2.*PHI50) ) / { 2.* ( PH I9 5-PH I 5 ) ) 

DO 608 L = 1 , 6 

IF(TLFS (L) .GE.FSK) GO TO 609 
6 C 8 CONTINUE 
609 IFSKTL=ITLFSC ( L ) 

Il=3*( IFSKTL-1 ) +1 
12 = 1 1+2 

WRITE (6 , 14 2) F SK, ( SKEW( I I ) , I 1 = 1 1 , 12) 

142 FORMAT (3 IX, ' SKEWNESS* , F8 .2,4X,3A8) 

C COMPUTE KURT CS I S AND DETERMINE CATEGORY 

FKG= ( PH I 9 5- PHI 5 ) / ( 2 . 44* ( PH 1 75-.PH 12 5 ) ) 

DO 612 L = 1 , 6 

IF (TLFK ( L) . GE. FKG) GO TO 613 

612 CONTINUE 

613 I F KT L= L- 1 

IF ( IFKTL.EQ. 0) I FKTL=1 
1 1 = 3* ( IFKTL-U+l 
12=11+2 

WRITE(6,162) F KG , ( KURT (ID ,11 = 11,12) 

162 FORMAT ( 3 IX, 'KURTOSI S’ , F 8 . 2 , 4X , 3A 8 ) 

1000 CONTINUE 



6h 



RETURN 

END 
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